On electrons and hydrogen-bond connectivity in liquid water 
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The network connectivity in liquid water is revised in terms of electronic signatures of hydrogen 
bonds (HBs) instead of geometric criteria, in view of recent X-ray absorption studies. The analysis 
is based on ab initio molecular-dynamics simulations at ambient conditions. Even if instantaneous 
thread-like structures are observed in the electronic network, they continuously reshape in oscilla- 
tions reminiscent of the r and t modes in ice (r^^lTO fs). However, two water molecules initially 
joint by a HB remain effectively bound over many periods regardless of its electronic signature. 



Water is an extremely intriguing liquid that contin- 
ues to excite the interest of scientists in many disci- 
plines. Many of its anomalous properties 0, origi- 
nate in the hydrogen bonds (HBs) among water molecules 
0, Q . The concept of a network liquid emerges naturally 
from this HB connectivity, an intuitive image that has 
provided interesting insights into the properties of wa- 
ter 0, Q . Direct structural information to characterize 
such network structure is hard to obtain experimentally. 
Diffraction techniques 0,0 offer radial distribution func- 
tions (RDFs) very naturally, but rely on reverse Monte 
Carlo techniques using force-field models to obtain fur- 
ther structural information |^ . Spectroscopic probes pro- 
vide a rich source of complementary information. X- ray 
emission Ull^ (XES) and X-ray absorption HI E 111 
(XAS) spectroscopies explore the electronic states of the 
liquid right below and above the Fermi level, respectively. 
In particular, the work by Wernet et al. [l^ has recently 
introduced an extremely interesting new component into 
the study of liquid water, by relating a pre-edge feature in 
the XAS spectra with broken HBs. The authors propose 
to determine connectivity by looking at an electronic- 
structure signature of the HBs. Their conclusion is dar- 
ing: the average coordination in liquid water would be 
^2 instead of the previously accepted value slightly un- 
der 4, displaying a filamentous picture, instead of the 
distorted, partly broken and fluctuating tetrahedral net- 
work described in so many papers before 0, 0, Q, • 

Is it true? This would be the wrong question to ask. 
The kind of network depends on the definition of the 
hydrogen bond, furthermore, on deciding whether two 
given molecules in a given configuration are bonded or 
not. There is no direct physical HB observable and there 
is arbitrariness in the choice of what is actually measured. 
Instead, we address the question of how relevant is the 
newly proposed network image for the description of the 
liquid in the sense of the insights it offers. The conven- 
tional criterion for HB is based on geometric consider- 
ations: an oxygen-oxygen distance within the first peak 
of the 0-0 RDF, and an upper critical bend angle a (see 
Fig. This "geometric" definition is based on total 



energy considerations in contrast to the newly proposed 
"electronic" one. In this paper we explore the adequacy 
of the newly proposed probe, including its time scale, vir- 
tually instantaneous as compared with atomic motions. 

Electronic structure calculations have been performed 
based on density- functional theory (DFT), within the 
generalized-gradient approximation (BLYP) 0, ll^ . 
The Siesta method is used 0, with a basis set of 
atomic orbitals at the double-^" polarized level [2^. For 
liquid water at ambient conditions, ab initio molecular- 
dynamics (AIMD) simulations have been performed in 
the microcanonical ensemble, based on the DFT forces 
and the Born-Oppenheimer approximation. Further de- 
tails are found in Ref. XAS spectra have been 
calculated for selected configurations (see below). The 
pseudo-atomic orbitals in the basis set have been PAW- 
transformed [i^ into all-electron atomic orbitals for cal- 
culating matrix elements. For our basis set, neglecting 
the very small inter-molecular matrix elements was found 
to give an adequate approximation for the purposes of 
this paper. The strong excitonic effect introduced by the 
attraction between the core hole and the excited electron 
is estimated in the Z -\- 1 approximation [22] • 

Notwithstanding the importance of the XAS experi- 
mental data for this and other purposes, the probe pro- 
vides a rather indirect measure of the electronics of the 
hydrogen bond, not least because of the mentioned exci- 
tonic effect. This and other difficulties [2^. (broaden- 
ing, alignments) make it very difficult to obtain quanti- 
tative comparisons for the liquid phase. We find it more 
useful for our purposes to use a ground-state probe that 
we validate against XAS data in cleaner systems. This 
validation is two-sided. On one hand we test our cho- 
sen probe, on the other, we test the extent to which the 
XAS probe reflects the properties of the electronic ground 
state. 

A very natural choice for describing electronic bonding 
within our method is bond order. In its simplest defini- 
tion [2^, the bond order between two atoms, 1 and 2, is 
Qi,2 = Y^Il Y^t P^^ '^1^1^^ where fi {v) sums over the basis 
functions associated to atom 1 (2), and p and S are the 
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density and overlap matrices, respectively. Bond orders 
depend on the choice of basis, and their arbitrariness has 
been described at length (see ^2^ and refe. therein). It is 
relative changes in that value what we use in this work, 
and these are shown below to be meaningful enough to 
support its conclusions. Mulliken's bond orders are com- 
pared with Mayer's jl^, from which the same conclu- 
sions are drawn. We believe that any other electronic- 
structure signature of the bonding [2^ would reflect the 
same physics. 

There has been a controversy m m s EH on the 

covalence of the HB and on its bonding or antibonding 
character, which would seem to affect our choice of elec- 
tronic probe. The electronic characteristics of the HB 
were nicely illustrated with maximally localized Wan- 
nier functions They can also be described in terms 
of an intramolecular polarization (rehybridization within 
the molecules) and intermolecular polarization or charge 
transfer (admixture of orbitals of different molecules), if 
using the language of single- molecule orbitals [23 ^ of 
atomic orbitals, in either mono- or multi-determinantal 
wave- functions [2^. The physics of the interactions be- 
hind the HB is, however, quite clear if one avoids the se- 
mantic problems that have been partly behind that con- 
troversy. In a typical HB there is an important electro- 
static attraction that dominates the energetics 
addition, there is a deformation of the electronic cloud 
around O's accepting lone pair towards the donating H, 
in response to the fleld generated by the latter. There 
are other effects (e.g. quantum fluctuations of the pro- 
tons) that are energetically less significant in principle, 
but could still be important for liquid water. They are be- 
yond the scope of this paper. The mentioned polarization 
of the lone pair happens at the expense of a slight con- 
traction of the electron cloud involved in the 0-H bond 
within the donor molecule, due to Pauli exclusion |3n| . 
This last deformation is the one behind the antibond- 
ing character of the HB discussed in the literature |2^ . 
However, the original deformation of the O's lone-pair to- 
wards the proton remains clearly bonding. In this study 
we thus concentrate on Qoh, between the donor H and 
the acceptor O. 

In order to compare XAS and Qoh we have calculated 
XAS spectra for a four-layer (001) ice slab, with the sur- 
face molecules exposing non-donated protons (dangling 
HBs). Repeated slabs are separated by 15 A of vacuum. 
An extra water molecule has been placed on top of one of 
the surface molecules, and its position has been varied. 
Figure shows the XAS calculated for the donor surface 
molecule and Qoh for that HB as a function of HB O- 
H distance for two different values of the flap angle, (j) 
(see Fig.^. (j) ^ 40° corresponds to an ideal tetrahedral 
arrangement as in ice. A zero flap angle gives a rather 
unfavorable situation, since the proton faces midway be- 
tween both lone pairs of O (in the nodal plane of the 
acceptor's HOMO orbital, of C2v's Bi character), which 
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FIG. 1: XAS peak intensity ratio /edge //pre (a), and bond 
order Qoh (b) versus intermolecular distance for two values 
of the flap, (j) = 40° and 0°, for the HB between an acceptor 
molecule and a clean (001) ice surface. 



partly inhibits the electronic deformation. 

We focus on two distinctive features of the XAS spec- 
trum directly related to the pre-edge observed in the ex- 
periments: (i) the relative intensity of the peaks related 
to edge and pre-edge features, and (ii) the energy differ- 
ence between both peaks AE. As discussed by Cavalleri 
et al. [2^, the spectral intensity is taken from the Z ap- 
proximation (initial-state), while AE is obtained from 
the Z -\- 1 approximation, given its sensitivity to flnal- 
state effects. Fig. ^ shows that = 0° produces a pre- 
edge twice as large as (/) = 40° for the same doo. This 
remarkable effect is closely reproduced by Qoh, as well as 
the distance dependence. AE increases with distance in 
a similar manner (not shown), also well replicated by the 
bond order (the effect of the flap angle is less noticeable 
in this case). A detailed study of the dependence of Qoh 
on intermolecular geometry in a water pair can be found 
in Ref. [SJ. It is important to note that both magnitudes 
(XAS and Qoh) agree in not displaying any obvious fea- 
ture (discontinuity, zero, minimum) that would deflne a 
natural threshold for HB breaking. We will thus refrain 
from establishing an arbitrary criterion for the moment 
and explore what can be learnt independent of it. 

The Qoh's are then calculated in an AIMD simulation 
of liquid water for all water pairs within a flrst coordina- 
tion shell, as deflned by the flrst peak in the 0-0 RDF. 
The flrst point that becomes apparent is that every wa- 
ter molecule is mainly donating one strong HB, while 
the second bond order is 2.2 times weaker in average, 
partly supportin g th e one-dimensional network picture 
proposed in Ref. [13. The asymmetry is, however, not 
extreme, as can be seen in the distribution of strong to 
weak bond-order ratio in Fig. [21 
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FIG. 2: (a) Distribution of strong to weak bond-order ratio 
in the liquid during a 25 ps AIMD run. (b) Evolution of the 
donating Qoh's for a given molecule [zoom of the inset]. 

Fig. I2tb) shows a fragment of the time evolution of 
the QoH on the two donated HBs for a representative 
molecule (the complete trajectory is shown in the inset). 
The figure shows a clear oscillation with a period of ~ 170 
fs. It would correspond to the intermolecular r and t vi- 
brations in ice (hindered rotation and translation) |3^ . 
An intermolecular 0-H under-damped oscillation with a 
period of ^170 fs has been indeed directly observed in liq- 
uid water using ultrafast infrared spectroscopy. The 
figure shows that the oscillation in Qqh is more pro- 
nounced than what expected from Fig. [2 (a), with lower 
values 5 times smaller than the higher ones in average. 
It shows that an anti-phase vibration of both donating 
HBs only accounts for part of the oscillation, the remain- 
ing part coming from the in-phase vibration that would 
weaken (strengthen) both donated HBs simultaneously. 
The figure also shows that a single HB survives many 
such extreme cycles before breaking (the average life-time 
for a HB is a fewps). 

HBs with low enough values of instantaneous Qoh are 
certainly contributing to the measured XAS pre-edge. 
The XAS probe is thus reflecting a very pronounced 
electronic effect, a swinging of the electron cloud of the 
lone pair following the "flapping" and stretching of HBs, 
which appears to XAS as if many HBs were broken. Con- 
sidering the anti-phase component, that behavior would 
give rise to a pulsating ID filament-like network image, 
relevant for electron dynamics. It is, however, irrelevant 
to the liquid dynamics, since it is apparent in Fig [2l that 
the molecules in a HB are still effectively bound even 
when the electron deformation is very small. The in- 
clusion of quantum fluctuations to the nuclear dynamics 
could even exaggerate the beating effect, since quantum 
and thermal fluctuations are comparable in scale (the 



zero point motion of a 170 fs oscillation is 12 meV). 

Figure ini shows the distribution of Qoh for all the wa- 
ter pairs in the liquid with roo < 3.5 A, for two differ- 
ent definitions of bond order. In both cases there is a 
clear minimum at small values of the bond order, with 
~25 % HBs below that threshold. It offers a natural crite- 
rion (albeit still arbitrary) for determining the presence 
of a HB, namely, Qoh > Qqh HI- 

Interestingly, the 

width of the distribution for the HBs (2AQoh = 0.030, 
Qou = 0.023) is comparable to that for a single oscillat- 
ing HB through time (the one in Fig. El 2AQoh = 0.026, 
OoH = 0.020). This means that the very different HB 
strengths are not so much due to different configurations 
as the liquid flows (on the time scale of several ps, the 
HB lifetime) but remarkably related to these 170 fs vi- 
brations. 
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FIG. 3: Distribution of bond orders in a 25 ps AIMD liquid 
simulation, using Mulliken's (a) and Mayer's (b). 

After monitoring the electronic deformation as signa- 
ture for connectivity, we finish this study by assessing 
the effect of such deformation on the liquid structure it- 
self. Having observed how the electronic cloud deforms 
in response to stretch and flap, it should be expected that 
the energetics would be affected by the flap, and so would 
the configurations visited in AIMD trajectories. In Fig.^ 
the distribution of flap angles obtained from AIMD and 
two classical non polarizable models are compared |3^ . 
The angles Oi and O2 (as defined in Fig. ^ within the 
first coordination shell are used to characterize both flap 
and twist. The distributions show clear differences. The 
main difference between SPCE and TIP5P is the fact 
that the latter puts negative charges around the center 
of the electron lone pairs, which induces a more realistic 
flap response. It is, however, exaggerated, since the po- 
larization is static. AIMD distributions reflect both the 
preference for the tetrahedral geometry and the flexibil- 
ity given by the dynamical response of the electron cloud. 
It has been argued 0, [l^, 113, that, in spite of these 
differences, AIMD lines up with any force field so far 
to produce qualitatively wrong configurational sampling 
of liquid water. It is certainly true that there are clear 
shortcomings in our scheme (BLYP approximations for 
exchange and correlation, neglect of protonic quantum 
effects) as well as in others. What we propose here, how- 
ever, is a highly plausible explanation of the experimental 
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results that does not imply the paradigm shift proposed 
by these authors. 




FIG. 4: Distribution of Oi and O2 angles (as in Fig.Q) in liquid 
water, for (a) SPCE ^ and (b) TIP5P 14] force fields, and 
(c) AIMD. Darker regions indicate higher values. The 61 = 62 
diagonal represents pure flap. Moving normal to that diagonal 
represents twist. The maximum value along the diagonal is 
for zero flap (abrupt cut for rigid molecules), (d) Isosurface 
of the AIMD distribution of donating H's around a molecule. 

In summary, using our electronic probe for HB connec- 
tivity we observe what could be described as ID filamen- 
tous structures, but they are pulsating in a 170 fs period, 
the geometric connectivity surviving intact for many pe- 
riods and thus many reshapings of the filaments. Even 
if this image of pulsating filaments is not relevant for 
the description of the liquid, it is likely that questions 
addressing its electronic structure could benefit from it. 
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